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Abstract: Land degradation has a major impact on environmental and socio-economic sustainability. 
Scientific methods are necessary to monitor the risk of land degradation. In this study, the environmental 
sensitive area index (ESAT) was utilized to assess land degradation sensitivity and convergence analysis in 
Korla, a typical oasis city in Xinjiang of China, which is located on the northeast border of the Tarim Basin. 
A total of 18 indicators depicting soil, climate, vegetation, and management qualities were used to illustrate 
spatial-temporal patterns of land degradation sensitivity from 1994 to 2018. We investigated the causes of 
spatial convergence and divergence based on the Ordinary Least Squares (OLS) and Geographically 
Weighted Regression (GWR) models. The results show that the branch of the Tianshan Mountains and 
oasis plain had a low sensitivity to land degradation, while the Tarim Basin had a high risk of land 
degradation. More than two-thirds of the study area can be categorized as "critical" sensitivity classes. The 
largest percentage (32.6%) of fragile classes was observed for 2006. There was no significant change in 
insensitive or low-sensitivity areas, which accounted for less than 0.4% of the entire observation period. 
The ESAI of the four time periods (1994—1998, 1998-2006, 2006-2010, and 2010—2018) formed a series 
of convergence patterns. The convergence patterns of 1994—1998 and 1998—2006 can be explained by the 
government's efforts to "Returning Farmland to Forests" and other governance projects. In 2006—2010, the 
construction of afforested work intensified, but industrial development and human activities affected the 
convergence pattern. The pattern of convergence in most regions between 2010 and 2018 can be attributed 
to the government's implementation of a series of key ecological protection projects, which led to a decrease 
in sensitivity to land degradation. The results of this study altogether suggest that the ESAI convergence 
analysis is an effective early warning method for land degradation sensitivity. 
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1 Introduction 


Land degradation, which is affected by extreme climate events, rapid population growth, and 
unsustainable land management practices, poses a global environmental threat. It is a major 
challenge facing the survival and sustainable development of individuals and communities around 
the world (Abu and Tumeizi, 2012; Gibbs and Salmon, 2015; Dubovyk, 2017; Mao et al., 2018; 
Horion et al., 2019). Forms of severe land degradation have occurred in 10%—20% of the world's 
drylands, the consequences of which are estimated to directly affect about 25010° people in the 
developing world (Reynolds et al., 2007). Xinjiang Uygur Autonomous Region has the widest 
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distribution and largest area of desert lands in China (Zhang et al., 2014; Ran and Li, 2015). Korla, 
a typical oasis city in Xinjiang, has undergone a swift and strong economic development 
accompanied by both urban expansion and industrial production. This development has made the 
area's land supply and carrying capacity highly complex (Ma and Gao, 2017). The interactions 
among the Tarim River, Taklimakan Desert, and other systems further make the ecological 
environment of Korla fragile and increase the risk of land degradation (Li, 2018). 

Land degradation is a complex and continuous dynamic process (Wang et al., 2014). Monitoring 
land degradation requires consideration of multiple factors. Based on remote sensing and GIS 
technologies, researchers used various indicators to identify the areas of land degradation. Key 
variables include spectral biophysical indicators, environmental indicators, and socio-economic 
indicators (Xu et al., 2009; Lamchin et al., 2016). Previous studies have generally lacked 
standardized methodology, however, or have centered on limited indicators (Weinzierl et al., 2016; 
Riva et al., 2017; Easdale et al., 2019). To comprehensively and straightforwardly describe the 
process of land degradation, it is possible to combine the quantitative and qualitative methods with 
specific indicators (Sommer et al., 2011; Salunkhe et al., 2018). 

The environmental sensitive area index (ESAI) approach was developed in the context of 
identifying and assessing potentially degraded Mediterranean regions (Kosmas et al., 1999). It is 
implemented by considering significant factors resulting in land degradation, namely, soil, climate, 
vegetation, and management indicators. These indicators can facilitate a deeper understanding of 
land degradation sensitivity. Further, the ESAI approach has advantages of simple, flexible, and 
rapidly implementation, and it can be adjusted to different-level scales as per the conditions of the 
study area (Kosmas et al., 2003). As a widely recognized approach, the ESAI has been successfully 
applied to field sites in Europe, Mediterranean Basin, and Central Asia (Contador et al., 2009; 
Salvati and Zitti, 2009a; Basso et al., 2012; Pravalie et al., 2017; Kolios et al., 2018; Jiang et al., 
2019). However, there have been few previous studies on the land degradation sensitivity in oasis 
cities such as Korla. 

The combination of land degradation sensitivity and convergence analysis is an extension of 
previous research. The concept of "convergence" originates in economists' measurements of the 
development gap among countries. It is widely used to in the study of economic and industrial 
development (Gugler and Pfaffermayr, 2004; Salvati and Zitti, 2009b; Manca et al., 2014), but 
rarely in the convergence analysis of environmental quality. In the context of environmental 
research, convergence analysis was firstly used to monitor the sensitivity of land degradation in 
Italy, which provided empirical evidence for spatial convergence at the country scale (Salvati and 
Zitti, 2008). In this study, convergence is indicative of a negative correlation between the 
differences in ESAI over a particular time interval and at the initial observation time. The ESAI 
with the lowest value is closer to the mean than other ESAI (Tombolini et al., 2016). Spatial 
convergence and divergence are applicable in assessments of regional land degradation and may 
provide a reference for governments in managing areas where degradation has occurred. Therefore, 
convergence analysis that includes environmental indicators may provide an evidence of the driving 
factors on land degradation. Results of convergence analysis can also serve as an early warning 
signal for land degradation sensitivity (Neumayer, 2001). 

In this study, the ESAI of Korla from 1994 to 2018 was calculated using time-series, multi-source 
spatial data to determine the sensitivity of land degradation across the region in different years. The 
rationality of these results was verified by field investigation. The convergence of the time-series 
ESAI was used to identify "hotspots" in need of attention and to analyze the effects of the 
implementation of land management measures. The novelty of this work is its focus on typical oasis 
areas, its validation of the feasibility of the ESAI approach, and its establishment of suitable 
methods for evaluating land degradation sensitivity in the areas significantly impacted by human 
activities and natural factors. Most importantly, a spatial-temporal analysis of ESAI was conducted 
to identify areas that require a close attention. 

The goals of this study are to (1) monitor the quality index and spatial-temporal pattern of ESAT; 
(2) analyze the spatial convergence of ESAI at time intervals on a grid scale; and (3) determine the 
underlying mechanisms of the convergence and divergence patterns of ESAI. The results presented 
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in this study may provide a scientific reference for optimizing land management and preventing 
land degradation in Korla, as well as other similar regions in the world. 


2 Materials and methods 


2.1 Study area 


Korla City (41°10'48"—42°21'36'N, 85°14’ 10"—86°34'21"E), located in the central part of Xinjiang 
Uygur Autonomous Region, China, is situated at the northeastern border of the Tarim Basin and 
crosses the southern reaches of the Tarim River alluvial fan (Fig. 1). It is situated in the Konqi River 
Delta, forming a vast and fan-shaped oasis plain. The northern part of Korla is a branch of the 
Tianshan Mountains that forms an alluvial fan in the front of the mountains. Korla has a total area 
of 7268 km?, with the length of 127 km and width of 105 km (Song et al., 2013). Korla has a 
temperate continental arid climate, with annual mean temperature of 11°C, mean annual 
precipitation of 59 mm, total annual solar irradiance time of 2990 h, and mean annual evaporation 
of 2778 mm. The favorable water, soil, light, and heat resources in Korla are very suitable for the 
growth of cash crops such as pear, cotton, and tomato (Liu et al., 2014). 

At the end of 2018, the total population of Korla was 4.73x10° and the regional gross domestic 
product (GDP) was 6.16x10!° CNY (Statistic Bureau of Koral, 2018). Korla has gradually become 
an important modern city in Xinjiang. However, land degradation caused by desertification, 
salinization, over-exploitation of cultivated land, and droughts significantly influence the 
ecological environment of Korla. These problems threaten the coordinated development of 
economy, society, and environment in the city and its surrounding regions (Cai et al., 2015). There 
is an urgent need to effectively monitor land degradation sensitivity in this area. 
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Fig. 1 Landsat TM imagery of the study area (Bands 5-4-3) 
2.2 Data and indicators 


The United Kingdom's Ministry of Agriculture, Fisheries, and Food launched the ESAI approach 
in 1987 to encourage farmers and landowners to adopt environmentally and friendly land 
Management measures (Wilson, 1996). Through field surveys and both general and local 
knowledge of the environmental processes at work, Kosmas et al. (1999) improved the ESAI 
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approach and applied it successfully in the Mediterranean Desertification and Land Use Project 
(MEDALUS). 

According to the ESAI approach, we selected 18 indicators in this study to take into account four 
quality domains: soil, climate, vegetation, and management. Each quality indicator was calculated 
by the geometric mean of the relevant sub-indicators. The values of each sub-indicator were split 
into several classes and several scores were assigned accordingly. The thresholds of the scores were 
determined with reference to extensive field work by the MEDALUS and as per the actual situation 
of the study area (Kosmas et al., 1999). We assigned sensitivity scores between 1 (the lowest 
sensitivity) and 2 (the highest sensitivity) to each class based on the importance of its role in land 
degradation sensitivity and its relationship to the onset of the degradation process (Kosmas et al., 
2014). 

To ensure a strong spatial consistency across different indicators, we resampled the spatial 
resolution of all sub-indicator data to 300 m using the nearest neighborhood method and on the same 
projection system. A time period encompassing 24 years was assessed by observing the land 
degradation sensitivity at five specific time points: 1994, 1998, 2006, 2010, and 2018. A field survey 
on local land situations was carried out in August 2018. Some sites representing land conditions of 
sparse vegetation were selected to further analyze the sensitivity of land degradation. 

2.2.1 Soil quality index (SQI) 

Soil can maintain biological productivity and environmental quality while supporting the health of 
plants and animals in various ecosystems (Salvati et al., 2011). According to the ESAI approach, 
we considered six parameters (texture, slope, surface albedo, dryness, salinity, and moisture) in the 
assessment of soil quality (Table 1). Soil texture data were downloaded from the Chinese Resource 
and Environment Data Cloud Platform (http://www.resdc.cn). The resource was compiled based on 
a soil type map with a spatial scale of 1:1,000,000. Soil profile data were obtained from the Second 
Soil Census (Tang, 1989), where soil texture was divided by the contents of sand, silt, and clay. 
The "Advanced Spaceborne Thermal Emission and Reflection Radiometer" (ASTER) Global 
Digital Elevation Model Version 2 with a spatial resolution of 30 m was used to retrieve the slope. 
The exact dates on which Landsat images acquired were 24 August 1994, 4 September 1998, 10 
September 2006, 20 August 2010, and 26 August 2018. These Landsat imagery scenes were chosen 
because they had almost no clouds and their acquisition dates were close together. All Landsat 
imagery scenes were radiometrically calibrated and atmospherically corrected. Surface albedo 
(Eq. 1), dryness (Eq. 2), salinity (Eq. 3), and moisture (Eq. 4) data were retrieved using Landsat 
images with a spatial resolution of 30 m. The band calculation formulas are as follows: 

Surface albedo = 0.356 x B+ 0.130 x R + 0.373 x NIR + 0.085 x SWIR, + 0.072 x SWIR, —0.0018, (1) 


2xSWIR, MIR Gg 
(SWIR, + R)-(NIR+B) _ SWIR,+NIR_ NIR+R_G+SWIR, 
(SWIR,+R)+(NIR+B) _2xSWIR, NR , G 
SI +IBI SWIR,+NIR NIR+R G+SWIR, 
Dryness = = p (2) 
2 2 
Salinity=/BxR, (3) 


Moisture = 0.0315 x B + 0.2021 x G + 0.3102 x R + 0.1594 x NIR — 0.6806 x SWIR, — 0.6109 x SWIR.,, (4) 


where B is the blue band; R is the red band; NIR is the near-infrared band; SWIR; is the short-wave 
infrared band 1 (1550-1750 nm); SWIR: is the short-wave infrared band 2 (2090-2350 nm); SI is 
the bare soil index; IBI is the index-based built-up index; and G is the green band. 

The SQI was calculated as follows (Kosmas et al., 1999): 


SQI=(Texture x Slope x Surface albedo x Soil dryness x Soil salinity x Soil moisture)”. (5) 
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Table 1 Calculation parameters of soil quality index (SQI) 


Indicator Parameter Class Description Score 
1 Loamy, sandy clay loam, and clay loam 1.0 
2 Sandy clay and silt loam 2 
Texture 
3 Silt, clay, and silt clay 1.6 
4 Sand 2.0 
1 <6 1.0 
Slope (%) 2 6-18 ie 
ope (% 

h 3 18-35 15 
4 >35 2.0 
1 Somber 1.0 
Surface albedo 2 Moderately bright L5 
3 Bright 2.0 
SQI 1 Low 1.0 
2 Moderate 13) 

Soil dryness 
3 High 1.6 
4 Very high 2.0 
1 None or slight 1.0 
2 Moderate 1.4 

Soil salinity 
3 Severe iL 
4 Very severe 2.0 
1 Very high 1.0 
2 High 1.3 

Soil moisture 
3 Moderate Ly 
4 Low 2.0 


2.2.2 Climate quality index (CQI) 

Severe climatic conditions such as low precipitation, gusty winds, and extreme temperatures can 
exacerbate land degradation. In this study, we assessed the CQI according to four parameters: 
precipitation, aridity index, aspect, and wind speed (Table 2). Precipitation and wind speed were taken 
from annual data sets of the National Meteorological Information Center (http://data.cma.cn/). The 


Table 2 Calculation parameters of climate quality index (CQD) 


Indicator Parameter Class Description Score 

1 >1000 1.0 
2 650-1000 1.3 

Precipitation (mm) 
3 280-650 1.6 
4 <280 2.0 
1 Humid (>0.65) 1.0 
2 Dry sub-humid (0.50-0.65) T2 
Aridity index (mm/mm) B Semi-arid (0.20-0.50) 1.5 
CQI 4 Arid (0.03-0.20) iy 
5 Hyper-arid (<0.03) 2.0 
1 N, NE, NW, E, and flat areas 1.0 

Aspect 

2 S, SE, SW, and W 2.0 
1 <2.0 1.0 
2 2.0-3.5 13 

Wind speed (m/s) 
3 3.5-4.5 1.6 
4 >4.5 2.0 
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Australian National University Spline (ANUSPLINE) meteorological interpolation model (Liu et 
al., 2008) was used to interpolate the precipitation and wind speed data from meteorological 
stations of varying altitudes and locations during 1994—2018, and the spatial resolution was 300 m. 
The aridity index (AI) with the spatial resolution of 1 km? was collected from the Consortium for 
Spatial Information (CSI). The data are available at http://www.cgiar-csi.org. The AI was calculated 
using the equation: AI=P/PET, where P is the annual precipitation (mm) and PET is the annual 
potential evapotranspiration (mm). Aspect is an important parameter that affects the angle and time 
at which sunlight strikes the soil surface, which can result in varying degrees of vegetation humidity 
(Kosmas et al., 1999). The ASTER Global Digital Elevation Model Version 2 
(https://search.earthdata.nasa.gov) was used to calculate the aspect. 
Specifically, the CQI was calculated as follows (Kosmas et al., 1999): 

CQI=(Precipitation x Aridity index x Aspect x Wind speed)". (6) 
2.2.3 Vegetation quality index (VQI) 
Vegetation regulates climate, maintains water and soil, reduces unwanted effects of wind and sand, 
and protects farmlands, especially in arid and semi-arid areas. The VQI was assessed in this study 
through four indicators: fire risk, erosion protection, drought resistance, and vegetation coverage 
(Table 3). The land use/cover map was provided by the European Space Agency the Climate 
Change Initiative-Land Cover (ESACCI-LC), a project that delivers consistent global land 
use/cover maps at 300 m spatial resolution on an annual basis. The data are available at 
http://maps.elie.ucl.ac.be/CCI/viewer/. The land use/cover map for the period of 1994-2018 was 
used to gather the detailed fire risk, erosion protection, and drought resistance information 
comprising the VQI. We then divided the results into different classes according to the land 
use/cover type. 


Table 3 Calculation parameters of vegetation quality index (VQI) 


Indicator Parameter Class Description Score 
1 Agricultural crops, water, wetland, bare land, and urban areas 1.0 
cman 2 Sparse vegetation T3 
Fire risk 
3 Shrubland and grassland 1.6 
4 Evergreen forest and conifer forest 2.0 
1 Evergreen forest, conifer forest, and urban areas 1.0 
2 Shrubland and grassland i183 
Erosion protection 3 Sparse vegetation 1.6 
4 Agricultural crops 1.8 
5 Water, wetland, and bare land 2.0 
VQI 1 Water, wetland, and urban areas 1.0 
2 Evergreen forest and conifer forest 12 
aaa 3 Shrubland and grassland 1.4 
4 Agricultural crops 1.7 
5 Sparse vegetation and bare land 2.0 
1 >75 1.0 
o ; 2 50-75 1.3 
S 25-50 15 
4 10-25 1.8 
5 <10 2.0 


The dimidiate pixel model was used to calculate the vegetation coverage in each period (Eq. 7) 


via the normalized difference vegetation index (NDVI) (Eq. 8) to determine the vegetation status. 
NDVI- NDVI 


FVC= sil 100%, (7) 
NDVI,,, — NDVI 


soil 
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-ARR © 
+R 
where FVC is the vegetation coverage (%); NDVIsoi is the NDVI value of the pure soil cover pixel; 
NDVlveg is the NDVI value of the pure vegetation cover pixel; NIR is the near-infrared band; and 
R is the red band. 

The VQI was calculated as follows (Kosmas et al., 1999): 


NDVI 


VQI = (Fire risk x Erosion protection x Drought resistance x Vegetation coverage)’. (9) 


2.2.4 Management quality index (MQI) 

Various anthropogenic environmental pressures and different forms of land use and management 
may affect land degradation sensitivity. In this study, the MQI was calculated using GDP, 
population density, agricultural intensity, and policy enforcement information (Table 4). GDP and 
population density data were obtained from the Chinese Resource and Environment Data Cloud 
Platform (http://www.resdc.cn). This data set uses the administrative area as a basic statistical unit 
to spread the population and GDP information onto a kilometer-level grid cell for spatialization. 
The original data can then be re-projected, re-sampled, and clipped as necessary. The agricultural 
intensity and policy enforcement information used here were originated from land use/cover maps 
provided by the ESACCI-LC with a spatial resolution of 300 m. According to local agricultural 
systems, the agricultural intensity depended on the type of plant varieties, the mode of irrigation, 
the degree of mechanization, and the other related factors at work in the given area. We classified 
policy enforcement across three categories according to the degree of land protection (Table 4). 

The MQI was calculated as follows (Kosmas et al., 1999): 


MQI=(GDP x Population density x Agricultural intensity x Policy enforcement)’. (10) 
Table 4 Calculation parameters of management quality index (MQI) 
Indicator Parameter Class Description Score 
1 <500 1.0 
2 500-1000 1.2 
GDP (gross output 3 1000-1500 1.4 
value/km°) 4 1500-2000 1.6 
5 2000—4000 1.8 
6 >4000 2.0 
1 <20 1.0 
2 20-60 1.2 
Population density 3 60-100 14 
mo (people/km?) 4 100-150 1.6 
5 150-200 1.8 
6 >200 2.0 
1 Evergreen forest, conifer forest, urban areas, water, wetland, bare 1.0 
land, and sparse vegetation 
Agricultural intensity 2 Grassland and shrubland i5 
S} Agricultural crops 2.0 
1 Urban areas 1.0 
Policy enforcement 2 Agricultural crops, water, wetland, grassland, and shrubland 1.5 
3 Evergreen forest, conifer forest, sparse vegetation, and bare land 2.0 


2.2.5 Environmental sensitive area index (ESAI) 


The ESAI is given by the geometric mean of four quality indices described above. The formula is 
described as follows (Kosmas et al., 1999): 


ESAI=(SQI x CQI x VQI x MQ)”. (11) 
According to the scoring system proposed by the MEDALUS, ESAI values range from 1 (the 
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lowest sensitivity) to 2 (the highest sensitivity). ESAI values can cover as diverse a range of land 
sensitivity to degradation as possible. Table 5 shows that ESAI values may be divided into eight 
classes reflecting the most commonly used classification thresholds (Kosmas et al., 1999; Basso et 
al., 2000; Salvati et al., 2011; Pravalie et al., 2017). 


Table 5 Classes and corresponding ranges of environmental sensitive area index (ESAI) 


Score range Class Description 

<1.17 N (non-affected) Insensitive to land degradation 
1.17-1.22 P (potential) Low sensitivity to land degradation 
1.23-1.26 F1 (fragile 1) Moderately sensitive to land degradation 
1.27-1.32 F2 (fragile 2) Moderately sensitive to land degradation 
1.33-1.37 F3 (fragile 3) Moderately sensitive to land degradation 
1.38-1.41 C1 (critical 1) Highly or very highly sensitive to land degradation 
1.42-1.53 C2 (critical 2) Highly or very highly sensitive to land degradation 

>1.53 C3 (critical 3) Highly or very highly sensitive to land degradation 


2.3 Methodology 


2.3.1 Spatial-temporal variation trend surface analysis 

Trend surface analysis not only reveals the spatial variation trends of data in different directions but 
also provides the optimal trend surface through function fitting. The trend surface reflects the overall 
variation of spatial objects in the spatial region (Li et al., 2012). In this study, trend surface analysis 
was used to project ESAI values onto planes in the east-west direction and the north-south direction 
of the study area. Function fitting was performed to obtain an optimal function curve that represents 
the trend of these projection points corresponding to the overall spatial-temporal variation trend of 
ESAI values. The final quadratic trend surface model has sufficient fitting accuracy. Its binary 
regression function is described as follows (Wang and Zuo, 2015): 

Z =a, +a,X+d,y+a,x° +a,xy+a,y +6, (12) 
where the spatial coordinates (x, y) of the sample points are independent variables; Z is the dependent 
variable; ao, a1, a2, a3, da, and as are the polynomial coefficients; and e is the error term. 

2.3.2 Convergence analysis 
"Convergence", in this context, refers to a negative correlation between the change rate of ESAI 
values and the ESAI value of the initial year during a particular time interval. Different regions 
have different initial ESAI values, and low-value regions may have faster growth rates. As time 
progresses, all regions will converge to the mean ESAI value, i.e., "reaching convergence". The 
average value of ESAI and annual average rate of different intervals were calculated as discussed 
in detail below. 

The following ordinary least squares (OLS) model (Kramer and Donninger, 1987; Tu and Xia, 
2008; Afendras and Markatou, 2016) was used to measure the overall spatial convergence of ESAI 
in Korla: 


k 
Y= Bo +X B; Xi tE, (13) 
j=1 


where y; represents the change rate of ESAI at the corresponding time interval; Jo is a constant 
term; Jj is a regression coefficient; i is the number of observations; x; (j=1, 2, ..., k) represents the 
ESAI value of the i™ spatial domain in the j" year; and e; is a random error term that conforms to a 
normal distribution. The coefficient J; is statistically significant (P<0.01) based on F-tests. A 
positive coefficient indicates a divergent ESAI trend, while a negative coefficient indicates a 
convergent ESAI trend (Salvati and Zitti, 2009b). 

The geographically weighted regression (GWR) model (Brunsdon et al., 1996) was used to 
identify the local spatial convergence of ESAI in Korla. A fixed weight Gaussian function was 
selected here as a spatial weight function and the AIC (Akaike Information Criterion) method was 


202010.00026v1 


chinaXiv 


ChinaXivA ERAT! 


JOURNAL OF ARID LAND 


used to determine the bandwidth of the weight function (Wang et al., 2005). The ESAI regression 
coefficient was obtained by local regression using data from adjacent observation points, and it 
varied with geographic location. The convergence analysis was estimated as follows (McMillen, 
2004; Tu and Xia, 2008; Ogneva et al., 2009; Liu et al., 2019): 


k 
Vi = By(u;v,) + $2, (u vi); +E, (14) 
j=l 


where y; represents the change rate of ESAI at the corresponding time interval; fo is a constant 
term; (ui, vi) is the geographic center coordinate of the i spatial domain; fo(ui, vi) is a constant 
term; f,(ui, vi) is the j regression coefficient of the i spatial domain; x; (j=1, 2, ..., K) represents 
the ESAI value of the i® spatial domain in the j} year; and ¢; is an error that obeys an independent 
normal distribution with a mean of zero. A positive ESAI regression coefficient indicates a 
divergent ESAI trend in the corresponding region, while a negative ESAI regression coefficient 
indicates a convergent ESAI trend in the corresponding region (Tombolini et al., 2016). 


3 Results 


3.1 Spatial assessment of ESAI from 1994 to 2018 


Figure 2 shows the spatial distribution of ESAI in Korla. Regions susceptible to land degradation 
were scattered throughout the study area. Generally, low ESAI values were concentrated in the 
branch of the Tianshan Mountains and the oasis plains, while high ESAI values were concentrated 
in the Tarim Basin, alluvial fan, and Tarim River alluvial plain, which represented the regions most 
susceptible to land degradation in Korla. The sensitivity of land degradation in desert areas was 
higher than that in mountains and oasis plains as expected. Anthropogenic activity had a little 
impact on grasslands in the branch of the Tianshan Mountains, so any slight changes may be a 
response to climatic factors. 


85°30'E 86°00'E 86°30'E 85°30'E 86°00'E 86°30'E 85°30'E 86°00'E 86°30'E 
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Fig. 2 Spatial distributions of ESAI (environmental sensitive area index) in 1994 (a), 1998 (b), 2006 (c), 2010 
(d), and 2018 (e), and spatial distribution of mean ESAI from 1994 to 2018 (f), as represented by different sensitivity 
classes of ESAI. The description of sensitivity classes is shown in Table 5. 
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Table 6 shows the area percentage of each ESAI sensitivity class. Across the observation period, 
most of the study area fell into critical classes. The area percentage of critical classes accounted for 
85.42% of the total area in 1994, but in 2006 it fell sharply to 66.98% before increasing to 76.33% 
in 2018. The fragile classes were mainly concentrated in the branch of the Tianshan Mountains and 
the oasis plains. In 1994, 14.39% of the total area belonged to these classes. This percentage 
remained above 20.00% from 1998 to 2018. Area percentage of non-affected and potential classes 
showed no significantly changes and accounted for less than 0.40% of the total area. In general, the 
sensitivity of land degradation in Korla declined over the study period. 


Table6 Area percentages of ESAI sensitivity classes in different years 
Area percentage (%) 


Class 

1994 1998 2006 2010 2018 

Non-affected 0.04 0.04 0.03 0.00 0.02 
Potential 0.15 0.33 0.39 0.03 0.31 
Fragile 1 0.41 1.50 2.35 0.40 i 2 
Fragile 2 4.91 9.02 13.89 6.97 7.38 
Fragile 3 9.07 14.98 16.36 13.83 14.74 
Critical 1 A T 15.42 15.14 12.00 13.25 
Critical 2 53.75 52.64 49.96 47.75 45.49 
Critical 3 18.90 6.07 1.88 19.02 17.59 


Overall, the sensitivity of land degradation in the west of Korla was much higher than that in the 
east (Fig. 3). These four land conditions of sparse vegetation of Tamarix ramosissima (Fig. 3a), 
cotton intensive agricultural areas (Fig. 3b), severely saline-alkali soils (Fig. 3c), and extensive 
desert areas (Fig. 3d) corresponded to the fragile 3, fragile 1, critical 2, and critical 3 classes, 
respectively. It can be seen that the ESAI accurately reflects this condition: high ESAI values are 
related to harsh environmental conditions. The land degradation differences among different 
locations in the field survey were consistent with the corresponding ESAI classification results. 

The trend line of land degradation sensitivity according to spatial statistical analysis in ArcGIS is 
shown in Figure 4. The blue line indicates the land degradation sensitivity trend in the east-west 
direction and the green line represents the trend in the north-south direction. The trend lines in the 
north-south direction formed convex curves in certain periods, indicating that the spatial distribution 
of ESAI was high in the middle and low on both sides of the study area. The spatial distribution of 
ESAI was the opposite in the east-west direction, i.e., low in the middle and high on both sides. The 
trend line of land degradation sensitivity was steeper in 2006 and 2010 (Figs. 4c and d), which 
indicates that land degradation sensitivity not only presents spatial differences of Korla but also 
changes in state and speed in different years. 


3.2 Convergence analysis of ESAI 


The convergence patterns were used to study the relationship between ESAI at the beginning of the 
time intervals assessed (1994, 1998, 2006, and 2010) and the difference over four respective time 
intervals (1994-1998, 1998-2006, 2006-2010, and 2010-2018). The coefficients of the four time 
periods obtained by the OLS model were —0.0095, —0.0102, —0.0114, and —0.0081, respectively, 
which indicates that there is a convergence trend in each period and the land conditions in Korla 
are stabilizing over time. 

The coefficients and local R? values of the ESAI spatial convergence analysis based on the GWR 
model at four time intervals are shown in Figure 5. Land degradation sensitivity in 96.93% of the 
study area showed convergence, and the convergence effect of the oasis plains was higher than that 
of the other regions from 1994 to 1998 (Figs. 5a1 and a2). During the period 1998-2006, regions 
with highly convergence were concentrated in the southeast with high local R? values, which 
suggests that ESAI tends to converge in the southeast during those years (Figs. 5b1 and b2). From 
2006 to 2010, the positive ESAI coefficients showing divergence were concentrated in the branch 
of the Tianshan Mountains and the oasis-desert ecotone with low local R? values, while regions 
with a higher convergence increased compared to the previous year and were mainly distributed in 


202010.00026v1 


chinaXiv 


ChinaXive FRAT 


JOURNAL OF ARID LAND 


the oasis plains with high local R? values (Figs. 5c1 and c2). Up to 18.95% of the regions were 
divergent in the oasis-desert ecotone, alluvial fan, and Tarim River alluvial plain; and the convergent 
regions were concentrated in the oasis plains with high local R? values for the period 2010-2018 
(Figs. 5d1 and d2). In general, land degradation sensitivity tended to converge between 1994 and 
2018. The ESAI converged or diverged faster in regions with high land quality than in regions with 
low land quality. 


85°30'E 86°00'E 86°30'E 
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Fig. 3 Spatial distribution of ESAI in 2018 (represented by different sensitivity classes; a) and the field survey 
photos of land degradation sensitivity areas (b, Tamarix ramosissima; c, cotton plantations; d, extensive desert areas; 
and e, saline-alkali soil) 


—— East-west trend line 
— North-south trend line 


Fig.4 Spatial trend surface analysis of land degradation sensitivity in 1994 (a), 1998 (b), 2006 (c), 2010 (d), and 
2018 (e) 
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Fig. 5 Spatial distributions of coefficients and local R? values of the ESAI spatial convergence analysis from the 
geographically weighted regression (GWR) model over the time intervals of 1994-1998 (a1 and a2), 1998-2006 
(b1 and b2), 2006-2010 (c1 and c2), and 2010-2018 (d1 and d2) 


4 Discussion 


In this study, we assessed a total of 18 indicators depicting soil, climate, vegetation, and 
management qualities according to the ESAI approach. The land degradation sensitivity and spatial 
convergence analysis of ESAI information were used to determine the potential land degradation 
in Korla from 1994 to 2018. The spatial convergence analysis reflects a series of complex 
interactions between ecological conditions and socio-economic backgrounds at different time 
intervals. Convergence patterns were observed in most regions over sequential time intervals. The 
convergence process gradually consolidates gaps between vulnerable regions. The oasis plains of 
Korla, which exhibited high land quality, were subject to relatively high anthropogenic 
environmental pressures (e.g., population growth, urban expansion, and crop intensification). The 
Taklimakan Desert in the western part of Korla showed a low-quality natural capital due to its harsh 
ecological environment and high sensitivity to land degradation (Yan et al., 2015). 

The convergence pattern of ESAI during the period 1994—1998 reflected the impact of macro- 
control measures such as city planning and industrial development. To strengthen the ecological 
environment, the local government of Korla implemented the "Returning Farmland to Forest" 
Program to encourage more farmers to plant fruit trees such as pears and red dates, and to stoke the 
enthusiasm of farmers by increasing subsidies for releasing their farmlands to forests (Song et al., 
2013). In addition, Korla has successively launched the "Three-North" Shelterbelt Project, resulting 
in remarkable achievements in forestry development (Abasi and Kasimu, 2015). 

During the period 1998-2006, the Chinese government implemented the "Western Development 
Strategy", which significantly accelerated the development of Korla. The population shifted from 
rural areas toward cities as urbanization grew rapidly, leading to a gradual increase in land bearing 
capacity (Emin and Simayi, 2014). The promotion of advanced agricultural irrigation technology 
has also strengthened the agriculture sector of Korla considerably. The ecosystem has been 
indirectly protected by these measures and the land use structure has been continuously adjusted. 
Thus, the areas of various major land use types have changed significantly (Figs. 2b and c). Areas 
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of cultivated land, garden, forest, and under-construction land showed overall upward trends while 
area of unused land decreased (Miao et al., 2012). The convergence pattern of ESAI during the 
observation period appears to be influenced by the plans of urban development and governance 
project. 

During the period 2006-2010, the construction of afforestation intensified, especially the 
construction of the Kongqi River Tourism Scenic Belt Project, which greatly improved the public 
green space of the city. As shown in Figure 5c1, a convergence phenomenon emerged nearby the 
Kongqi River Basin. The rapidly developing industry has brought a certain load to the environment. 
Building sites increasingly occupied large areas of land and have reduced the available green space. 
The divergence pattern of the partial regions is shown in Figure 5c1, which reflects the impact of 
anthropogenic activity (Aizezi et al., 2014). 

From 2010 to 2018, both land urbanization and population growth resulted in an increase in the 
development of Korla (Ma and Gao, 2017). The divergence pattern appeared to be significantly 
affected by the scale of anthropogenic activities. The government adopted effective farmland 
protection measures and scientific urban planning layouts, increased land development and 
rectification efforts, and realized a balance of cultivated land occupation and compensation 
(Maimaiti et al., 2017). The city structure was continuously optimized and adjusted alongside 
breakthroughs in industrialization. Advanced urbanization has increased public awareness of 
cultivated land protection. In 2018, the government continued to implement key ecological projects 
such as the northern ecological protection corridor, natural forest protection, sand control, and soil 
and water conservation efforts to improve the protection of natural Populus euphratica forests 
(regarded as an ecological security barrier system) on the lower reaches of the Kongqi River, which 
has reduced the land degradation sensitivity of Korla. 


5 Conclusions 


In this study, the ESAI approach was used to assess the land degradation sensitivity in Korla from 
1994 to 2018. A convergence analysis on the ESAI was performed at different time intervals. The 
land degradation sensitivity of Korla not only presents spatial differences but also changes in state 
and speed in different years. In general, land degradation sensitivity in Korla declined over the 
years examined in this study. The ESAI values in four time periods showed a convergence pattern, 
suggesting that land conditions are growing increasingly stable in Korla. The ESAI converged or 
diverged faster in areas with high land quality than those with low land quality. 

Limited data availability and spatial resolution restricted the analyses presented in this paper in 
regards to the sensitivity of land degradation at five time nodes. In the future, under permission for 
further data acquisition, it will be possible to analyze longer-term series with finer temporal 
resolution. The precise mechanism of land degradation also merits further research to elucidate the 
impact of various indicators on land degradation sensitivity. It is yet necessary to secure a workable 
reference for further research on the dynamic changes and convergence of land degradation 
sensitivity in arid areas. 
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